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The literature of marine invertebrate zoology contains relatively few references 
to genetic studies. Buzzati-Traverso, in a summary of accomplishments and a 
prospective future of marine biology, declared that knowledge of the genetics of 
marine organisms was in its infancy and that an understanding of evolutionary 
processes in the ocean required much more knowledge in that area (Buzzati- 
Traverso, 1960). Only a few investigations have detected unit genes in the 
Mendelian sense in which alleles with contrasting phenotypic expression were 
shown to segregate in the progeny. A partial list of works of this type includes 
those on the amphipod Gammarus chevreuxi (Sexton and Clark, 1936), the 
copepod Tisbe reticulata (Bocquet, 1951; Battaglia, 1958), the isopod Jaera marina 
(Bocquet, 1953), the polychaete Pomotoceros triqueter (Feyn and Gjøen, 1954), 
and the ascidian Botryllus schlosseri (Sabaddin, 1959; Milkman, 1967). Most 
investigations began by noting the occurrence of color polymorphism in the species 
concerned and led to the elucidation of a simple genetic controlling mechanism. A 
few other aspects of marine genetics have received brief attention. The cytogenetics 
of the gastropod Thais lapillus in different environments was studied (Staiger, 
1954+), and the polygenic basis of sex ratios in the copepod Tigriopus was investi- 
gated (Ar-Rushdi, 1958). Barigozzi (1960) speculated on the nature of genetic 
systems in sessile marine organisms. Useful reviews of the state of marine genetics 
have been published (Buzzatti-Traverso, 1960; Ray, 1960; Battaglia, 1965). The 
chief conclusion to be gained is that genetic studies on marine animals have been 
extremely limited in scope and theoretical significance. 

Other than chromosome counts of a few species (Makino, 1951), nothing is 
known of the genetics of the Phylum Ectoprocta. Marine ectoprocts are abundant, 
sessile benthic animals that feed on phytoplankton by means of a ciliated crown 
of tentacles, the lophophore. The initial phase is a solitary planktonic larva, which 
soon settles. A colony is formed by successive asexual budding of new individuals 
from previously formed individuals. Thus each individual (zooid) of an ectoproct 
colony 1s genetically identical. 

Several fundamental aspects of the life cycle of species of ectoprocts are poorly 
known. A point of special confusion has involved the mode of fertilization because 
the majority of ectoprocts are hermaphroditic. According to Hyman (1959), 
most specialists believe that hermaphroditic ectoprocts are self-fertilizing. Silén 
(1966) and Bullivant (1967) concluded from observations of sperm release into 
sea water that at least some species may practice outbreeding. Resolution of the 
question of either autogamy or intra-colony selfing versus outbreeding, and the 
understanding of genetic systems in general have been impeded by two practical 
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difficulties. (1) Morphological variants of potential genetic origin have not been 
recognized in any ectoproct species: and (2) ectoprocts have not been successfully 
crossed in the laboratory. One or both of these difficulties apply equally well to 
many marie animals. 

For the last decade, and especially since 1966 when the papers of Hubby and 
Lewontin, and Lewontin and Hubby appeared, electrophoretic and histochemical- 
staining techniques have been utilized for the detection of genetic variants at single 
loci. Crude protein extracts are prepared from individual organisms and spotted 
onto starch or polyacrylamide gels. Polypeptides coded by each allele of a locus 
migrate a characteristic distance in an electric field applied across the gel. The 
mobility of the polypeptide depends on its electrical charge, size and contormation. 
Mutations that alter these properties of a polypeptide will alter its mobility; 
therefore the number and kinds of alleles at a locus can often be determined by the 
electrophoretic pattern displayed by the corresponding polypeptides after fixation 
and staining by the proper histochemical means. 

A considerable body of literature on electrophoretic variants of enzymes now 
exists and much of it concerns genetic variation (reviewed by Shaw, 19653). The 
use of electrophoresis in population genetics has been strikingly exploited in studies 
of Drosophila (for example, Hubby and Lewontin, 1966; Lewontin and Hubby, 
1966; Johnson, Kanapi. Richardson, Wheeler and Stone, 1966; Yarbrough and 
Kojima, 1967; Kojima and Yarbrough, 1967; Stone, Wheeler, Johnson and 
Kojima, 1968; O’Brien and MacIntyre, 1969). A beginning has been made in 
the application of these techniques to other organisms, most notably by Manwell 
in taxonomic problems of marine invertebrates (Manwell, 1966; Manwell and 
Baker, 1963, 1968; Manwell, Baker, Ashton and Corner, 1967). Much of the 
protem variation in his research was ascribed to genetic causes, but gene loci were 
not usually identified. 

We used electrophoretic techniques to identify gene loci m local populations 
of several species of marine ectoprocta. Results for the two species best known 
are reported here. Where our collections were sufficient, gene and genotype 
frequencies were determined for individual loci, and the amount of genetic variation 
was calculated. 

Bugula stolonifera is erect and foltaceous in growth form, and is lightly calcified, 
Whereas Schizoporella unicornis grows as an encrusting, highly calcified sheet 
(Schopf and Manheim, 1967, 1968). Each of these species is widely reported from 
the Atlantic Coast of North America, England, and Western and Southern Europe 
(Ryland, 1960, 1965; Maturo., 1966). 


NLATERIALS AND METHODS 


Iectoproct colonies were collected in the vicinity of Woods Hole, Massachusetts. 
Schisoporella unicornis was obtained from pilings of the Marine Biological Labora- 
tory (AIBL) dock facing into Vineyard Sound, from a floating dock in Green 
Pond, Massachusetts, 10 km east of Woods Hole, and from a dredge haul from 
Vineyard Sound about 35 km southwest of Woods Hole (see Fig. 6). Most of the 
Bugula stolonifera came from metal drums used as floats for the MBL supply 
department area in Eel Pond, Woods Hole. Supplementary collections were 
taken from Green Pond. 
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Colonies were kept alive in running seawater tanks. They were prepared 
for electrophoresis by first removing any tmtimately associated fauna, such as 

phipods and nematodes from Bugula, and overgrown barnacles, tunicates and 
bivales from Schisoporella. Portions of Bugula colonies containing 1000-1500 
living zooids were triturated in a tissue grinder in cold 123% sucrose buffer. 
\bout 500 zooids per colony of the much more heavily calcified Schizoporella were 
ground in sucrose-buffer in an Albers micromortar (Thomas Co.). Ground 
samples were centrifuged at room temperature for 15 minutes at 3500 rpm in 
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Figure 1. Diagrammatic sketch showing patterns of bands illustrating genetic loci: 
A. Simple diallelic system with an upper band homozygote, a heterozygote, and a lower band 
homozygote; B. Diallelic system with an upper band homozygote, a heterozygote with a 
hybrid band, and a lower band homozygote; C. Inhibition of middle zone without affecting upper 
or lower zones; D. Monomorphic middle zone which does not vary in position relative to zones 
of defined loci above and below it; E. Polymorphic middle zone which varies in position relative 
to zones of defined loci above and below it; F. Zone of several bands of possibly more than one 
locus but which cannot be further resolved by techniques at hand. Arrows denote the zone 
under consideration. 
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1 ml polyethylene centrifuge tubes and 10 pl of supernatant were pipetted onto 
the gels. 

Electrophoresis was performed with a vertical electrophoresis cell (E-C 
Corporation, Philadelphia, Pennsylvania) similar to the method of Hubby (1963) 
and Hubby and Lewontin (1966). However, a 63% polyacrylamide gel was 
used throughout because of its superior resolution of bands of malate dehydrogenase. 
Electrophoresis was conducted for 3 hours at 350-400 V, and 50-100 mA for most 
analyses. Recirculating water from an ice bath kept the buffer temperature near 
15° C. Usually 23 colonies were electrophoresed simultaneously. The 24th slot on 
the gel was reserved for Bovine serum albumin (Nutritional Biochemicals Corp.). 
The major, fast moving band of Bovine albumin served as a reference standard 
for distance in all analyses. 

Electrophoresis of nonspecific esterases was conducted using 0.2 ./ Tris-glycine 
buffer, pH 8.3. Gels were preincubated in 150 ml of 0.5 AJ boric acid at 5° C for 
one hour and then transferred to 150 ml of 0.1 A/ phosphate buffer solution, pH 
6.5 at room temperature, containing 40 mg of a-naphthyl acetate as substrate and 
75 mg of Fast Red TR Salt as dye coupler. Staining was complete in two hours, 
after which gels were photographed and described, and then preserved in 4% 
acetic acid. 

Esterase inhibition studies were conducted on 3—5 colonies with eserine sulfate, 
EDTA, and phenylmethylsulfonylfluoride (PMSF). In alternate substrate studies, 
40 ing B-naphthyl acetate and 0.3 ml e-naphthyl butyrate were substituted for the 
usual substrate. 

Determinations of malate dehydrogenase were made using 0.2 JJ Tris-glycine 
buffer. pH 8.3. Tris-borate buffer, pH 8.9, and Tris-EDTA-borate bufter, pH 38.9; 
were also tried but did not yield good resolution. After electrophoresis, gels were 
incubated in 150 ml of 0.003 JY Tris-HCl-buffer, pH 8.5 at room temperature, 
containing 20 mg of NAD, 60 mg of Nitro BT chloride and 3 mg of phenazine 
methosulfate. Bands were developed after 14-24 hours of incubation in the dark. 

Electrophoretic determinations of “leucine” amimopeptidase were made with 
0.1 Af Tris-borate buffer, pH 8.9, containing ImJ/ MgCle. Following a preincuba- 
tion period of 30 minutes in 125 1m] 0.1 M Tris-maleate buffer, pH 5.3, and 40 mg of 
Leucyl-8-naphthylamide at room temperature, 75 mg of Fast Black K salt (Pfaltz 
and Bauer, Co., New York) were added. Bands developed during the following 
1-2 hours. Occasionally it was necessary to remove the old staining solution and 
add fresh substrate and dve coupler to permit development of faint bands. 


IDENTIFICATION OF GENETIC POCI 


In order to obtain genetic information froin electrophoresis gels, a means must 
be found to estimate the number of gene loci and the number of alleles per locus 
from the pattern of enzymatically active bands. Much of the difficulty of interpre- 
tation is obviated if organisms suspected of carrying different allelic combinations 
at a locus are crossed under controlled laboratory conditions. The progeny are then 
scored for patterns and these are compared to expected genotypes. Most discoveries 
of electrophoretic variants under genetic control have been made in maize, 
Drosophila, Mus, and other organisms that are readily crossed (Shaw, 1965). 
Many marine invertebrates, and ectoprocts in particular, have never been success- 
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fully crossed in the laboratory. Therefore three criteria that do not involve 
crosses were utilized in this study in estimating the number of genetic loci (Fig. 1). 


(1) A codominant diallelic or triallelic locus in which each allele produces a 
single enzyme band of characteristic mobility is easily recognized (Fig. 1A). 
Heterozygotes have the bands of both alleles. Each heterozygote band ts usually 
lighter than the corresponding band of the homozygotes since it 1s produced by a 
single dose of the gene. This type of band pattern has been noted repeatedly in 
reference in Drosophila and other organisms (see, for example, Wright, 1963: 
Beckman and Johnson, 1964a; 1964b; Scandalios, 1969). Tt also occurs m the 
ectoproct species under study, particularly in esterase gels of Bugula stolonifera 
and LAP gels of both species. We regard such a pattern as a reliable criterion 
for the existence of a polymorphic genetic locus. 

In multimeric enzymes, polypeptide subunits may form one or more inter- 
mediate bands approximately evenly spaced between the bands of homozygote 
mobilities. Commonly a single intermediate band is found in heterozygotes (Fig. 
1B) which presents no problem in interpretation. However, the genetic relation- 
ships of multimeric enzymes manufactured by several mobility alleles can be 
complicated (for example, see Burns and Johnson, 1967). Complex loci of this 
type occur in the species Bugia turrita and B. simplex (unpublished data) but not 
in B. stolonifera and S. unicornis. 

(2) The use of specific inhibitors and alternate substrates also permits the 
recognition of loci. The rationale is that a band or a band complex differentially 
affected by these agents in comparison to adjacent zones probably includes poly- 
peptides that originate from a single locus (Fig. 1C). This conclusion does not 
extend to groups of bands with widely different mobilities that happen to respond 
to the same inhibitors or alternate substrates. 

Gels stained for nonspecific esterases with a-naphthyl acetate usually show 
numerous zones of activity that represent several enzymes. Alternate substrates 
and various inhibitors with partial specificity have been used to distinguish classes 
of esterases (Allen, 1961; Ogita and) Kasai, 1965; and “‘Mianwell ci alin 19678 
Eserine sulfate inhibits certain cholinesterases, and EDTA inhibits some aryl- 
esterases (Augustinsson, 1961). PMSF inhibits acetycholinesterase, trypsin and 
chymotrypsin (Gold, 1967). The alternate substrates -naphthyl acetate and 
a-naphthyl butyrate are preferentially acted upon by certain esterases. 

(3) Bands are not randomly scattered over the gels, but are usually localized 
into discrete zones. A zone is considered representative of a single locus if variation 
within it appears independent of other zones of the gel (Fig. 1D). Theretore, a 
zone of activity appearing as a band of uniform mobility is regarded as mono- 
morphic if it is adjacent to zones definitely identified as polymorphic (Hubby and 
Lewontin, 1966). Similarly, zones containing bands of variable mobility but 
difficult to interpret genetically are considered polymorphic if they are adjacent to 
zones identified as either monomorphic or polymorphic by other criteria (Fig. 1E). 

We have no criterion for distinguishing loci where series of closely adjacent 
bands of uniform mobility are present other than the use of inhibitors and alternate 
substrates as in the case of esterases (Fig. 1F). A series of bands may represent 
several independent monomorphic loci, a sharing of subunits in a multimeric 
system or an artifact of conditions of the analysis. In this study, series of adjacent 
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bands of uniform mobility have been conservatively regarded as belonging to the 
Same locis, 

Based on these determinations we have employed a terminology for enzymes 
and genetic loci based principally on the system of Hubby and Lewontin (1966). 
For example, in B. stolonifera there are 6 zones of esterase activity under the 
control of separate genetic loci (Fig. 2). These enzymes are designated esterase-1 
through esterase-6 in order of increasing mobility and distance from the origin. 
The corresponding loci are termed E-1 through E-6. Thus, the fourth zone of 
esterase activity from the origin is the enzyme esterase-+ and it is controlled by 
the locus E-4. The esterase-4 enzyme has two mobility variants that migrate 
distances of 1.01 and 1.04 relative to that of the bovine albumin standard. The 
mobility variants belong to two alleles of the E-4 locus, and are designated E-+'° 
and E-t: (Table T). i 

Alleles of some loci produce two or more enzyme bands when homozygous. 
The mobility of such alleles is designated by a single superscript which is the mid- 
point between the outermost bands. Thus, the monomorphic locus M-2:*° of Bugula 
stolonifera corresponds to three bands of malate dehydrogenase ranging between 
0.75 and 0.95 in mobility. 

The number of alleles of the polymorphic loci E-3 in Schisoporella unicornis 
and M-1 in B. stolonifera is as yet unknown. Superscripts have not been applied 
to these loci. 


RESULTS 
Esterases 
Sehisoporella unicornis and Bugula stolonifera are quite unlike in the number 


and mobilities of zones of esterase activity (Fig. 2, Table I). The simplest band 
pattern occurs in S. wnicornis, in which there are 4 major zones of activity, 


TABLE | 


Locus and allele nomenclature for esterases and biochemical characterization 
of loct in Schizoporella unicornis and Bugula stolonifera 
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Ficure 3. Plot of areal distribution of genotypes of Bugula stolonifera (locus E-5) on a 
floating drum at Eel Pond. This pattern of random distribution of a genotype is typical for 
other loci on this drug and on others. 


corresponding to + gene loci, 3 monomorphic and 1 polymorphic and diallelic 
(Table I). Esterase-1, the product of the monomorphic E-1 locus, is inhibited by 
1 x 10+- Af eserine sulfate, aud does not appear when fB-naphthyl acetate and 
a-naphthyl butyrate are the substrates. Enzymes of the other loci are unaffected by 
the inhibitors employed and show little differential activity on alternate substrates. 

Esterase band patterns are more complex in B. stolonitfera. The 6 major zones 
of esterase activity are each controlled by a separate locus, E-1] through E-6. Four 
of the loci are polymorphic and diallelic (Table 1). As in S. unicornis, esterase-1 
is inhibited by eserine sulfate and does not utilize B-naphthyl acetate and a-naphthy! 
butyrate. Esterase 3 is inhibited by PMSF and the two bands of esterase-6 are 
greatly clarified by EDTA. 


Malate dehydrogenase 


The staining procedure for malate dehydrogenase (MDH) has greater 
specificity than that for esterases. Nevertheless, several MDH enzymes with 
different coenzyme and divalent cation requirements are known (Moore and 
Villee, 1963; Manwell ef al., 1967). In our study, attempts to detect NADP- 
and Mg*t-dependent MDH were unsuccessful. NAD-dependent MDH., however, 
was detected in both species. 
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In Schizoporella unicornis, MDH activity appears in a single band of uniform 
mobility (Fig. 2). It is the product of the monomorphic locus M-1. 

The fastest migrating three bands of Bugula stolonifera have constant mobilities 
in all colonies examined. We conservatively regard them as controlled by the 
same monomorphic locus, M-2. [n contrast, the bands of lesser mobility occur 
as two variants. In vanant A (Fig. 2, A, B), two sharp bands with mobilities of 
0.58 and’ 0:68 are present.” In variant 6b (Fie. 2, C) the two bands om Miani A 
are fainter and are joined by three additional, evenly spaced bands. This 5-banded 
pattern of variant A resembles the band pattern of the heterozygote of a 
tetramer. It is possible, however, that the variation is non-genetic, as has been 
noted for various dehydrogenases (see Jacobson, 1968; Thornber, Oliver and 
Scutt, 1968; and O'Brien and MacIntyre, 1969). On the other hand, the two 
variants remain clearly recognizable m spite of widely different amounts of crude 


TABLE [Í 


Locus and allele nomenclature for malate dehydrogenase and “leucine” 
aminopeptides in Schizoporella unicornis and Bugula stolonifera 











DEN Locus Monomorphic or Nomenclature of 
; S S z 
Type of protein vets Polymorphic locus and alleles 
Malate dehydrogenase Schisoporella unicornis Monomorphic M-1-9 
Bugula stolonifera Polymorphic M-1 
Monomorphic M-2-%5 
“Leucine” aminopeptidase Schizoporella unicornis Monomorphic Lap-1-*! 
Monomorphic Lap-2 -$$ 
Polymorphic Lapo 
Lap-3-%8 
Bugula stolonifera Monomorphic Lap-1 -$ 
Monomorphic Lap 2 
Polymorphic Lap-3-3 
Lopo 
Lap-3% 














protein introduced into the gel slots, different parts of the colony electrophoresed 
(growing apical zooids, reproducing zooids, and degenerating zooids (brown body 
regions) ), and differences in pH and ionic strength of the buffer. This is good 
evidence that they have a genetic basis. If the variants are genetic, probable 
differences in gene frequency exist between Eel Pond and Green Pond. Only 4 
of 28 colonies (or 14%) from Eel Pond are variant B whereas 11 of 26 (42%) 
from Green Pond are variant b. Variants A and B of B. stolonifera are treated 
as representative of one polymorphic locus, M-I. 


“Leucine” aminopeptidase 


Zones of activity for “leucine” aminopeptidase (LAP) are similar on gels of 
Schizoporella unicornis and Bugula stolonifera. Both species possess three loci 
controlling LAP, two monomorphic and one polymorphic (Table II). The 
enzymes of the Lap-1 loci of both species are similar in mobility (Lap-1-** in 
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Jaane I 


Amount of polymorphism in esterases, malate dehydrogenase, and “leucine” 
aminopeptidase for Schizoporella unicornis and Bugula stolonifera 


pe | 











Schizoporella unicornis | Bugula stolonifera Totals |Per cent polymorphism 

Esterases | | 

Monomorphic S 2 5 

Polymorphic 1 4 5 50.0. 
Malate dehydrogenase | 

Monomorphic 1 1 2 | 

Polymorphic 0 1 | Lay 33S 
“Leucine” aminopeptidase | 

NIonomorphic 2 2 4 

Polymorphic 1 | R 33 
Totals 02 56 11 8 42.1 
Per cent polymorphism 29.0 54.5 | 42.1 








S. unicornis and Lap-1?5 m B. stolonifera). Lap-2 enzymes are quite different 
11 mobilities and staining properties. The Lap-3 locus in S. wsticornis has two 
alleles and the corresponding locus of B. stolonifera has three. Once again, the 
protein mobilities are similar (mean mobility of 0.96 in S. wnieornis and 0.88 in 
B. stolonifera). The similarity in banding pattern of Lap-l and Lap-3 may 
indicate the presence of homologous proteins that are conservative to changes in 
electrophoretic mobility. 


Gene and genotype frequencies and polymor phisi 


Genotype frequencies have been determined for 5 separate polymorphic loci 
(Table IV). Data for Lap-3 of Sciizoporella wincormis are given for 3 localities 
making a total of 7 determmations. Gene frequencies were calculated from geno- 
type frequencies by the formulas: p= P + 4H, and q= 0+ 4H, where p and q 


TABLE IV 


Summary of gene and genotype frequencies and chi-square agreement with Hardy- 
Weinberg equiktbrium values for Schizoporella unicornis and Bugula stolonifera 














eeu Genotype frequency 
Tp es: Sample 7 a i ae. oe e 
Species ocality k ocus A p of genotype frequency 
ri j Zii Ea L High with Hardy-Weinberg 
mo- | mo- |mobility | Hetero- mobility 
bility | bility | homo- | zygotes | homo- 
allele | allele | zygotes | zygotes 
Schizoporella | MBL Dock | 50 Lap-3 0.720] 0.280 25 22 3 no significant 
24922COF RLS | difference 
Green Pond 43 Lap-3 0.756 | 0.244 25 15 3 no significant 
difference 
Vineyard Sound 33 Lap-3 0.485 | 0.515 9 14 | 10 no significant 
difference 
Bugula Eel Pond 44 E-2 0.864 | 0.136 34 8 i 2 no significant 
stolonifera | difference 
Eel Pond 43 E-3 0.602 | 0.398 20 12 11 P < 0.01 
Eel Pond 44 E-4 0.955 | 0.045 40 4 0 no significant 
difference 
Eel Pond 44 E-5 0.693 | 0.307 20 21 3 no significant 
| difference 
f 
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Fıcure 4. Distribution of genotypes of “leucine” aminopeptidase (locus Lap-3) in colonies 
of Schisoporella unicornis occurring on outer floating dock of boat moorings at Green Pond, 
Massachusetts. White area represents the surface available for colonization. Note that the 
boat mooring side has a greater frequency of heterozygotes than the open water side. The 
difference is significant at the 5 per cent level. 
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are gene frequencies and P, H and Q are the trequencies of homozygote, heterozy- 
gote and other homozygote classes respectively (Falconer, 1960). The Hardy- 
Weinberg equilibrium distribution of genotypes was calculated from gene fre- 
quencies for the 5 loci. Genotypes obtained were compared to expected values 
for goodness-of-fit (Table IV). Only the genotype frequency of the E-3 locus 
in Bugula stolonifera differs significantly from Hardy-Weinberg equilibrium values. 
At this locus there is a significant deficiency of heterozygotes (P < 0.01). 

In Schizoporella unicornis 2 of 8 or 25.0% of loci are polymorphic and in 
Bugula stolonifera 6 of 11 or 54.5% of loci are polymorphic. Mean polymorphism 
is 42.1% (Table IIT). Combined data from both species show 5 esterase loci 
COD. IMIDE. locus (33.39%. and 2 C Aci (93.30) to Dapon nigepmc 


Geographic variation in gene and genotype frequencies 


Local populations or demes of Bugula stolonifera in Eel Pond extend over 
areas of at least 3.6 m°. Genotypes of the E-5 locus are randomly distributed over 
Piisanen ( Fig. 3). 

Collections were made from the floating dock at Green Pond and the pilings at 
the MBL dock which show the distribution of Lap-3 genotypes of Schisoporella 
unicornis over areas of 11-18 m? (Figs. 4,5). At Green Pond (Fig. 4) there is a 
possible bias of heterozygotes (x° = 4.26, P < 0.05) on the side of the dock which 
is part of a motor boat mooring. The water on this side is more polluted than 
that of the free harbor. In contrast, the distribution of genotypes in a sample 
collected from several pilings of the MBL dock appears to be random (Fig. 5). 
The depth distribution on the pilings, however, does suggest a non-random dis- 
tribution of genotypes (Fig. 5). The homozygote Lap-3-9'/-%* is absent in 22 
near surface samples, occurs in 1 of 14 samples of intermediate depth, but accounts 
for 2 of the 4 deepest samples. 

A third geographic comparison is available for populations separated on the 
order of kilometers (Fig. 6). S. unicornis from Vineyard Sound dredgings differs 
significantly from the MBL dock and Green Pond populations in gene frequency 
aeethe Lap-3 locus (x° = 9.3, P < 0.01, in a comparison of themWBl dock and 
Vineyard Sound populations). The frequency of Lap 3° averages about 24% 
higher in the Vineyard Sound collection. Samples of MBL dock and Green Pond 
populations do not show significant differences. 


DISCUSSION 


The most important conclusion we can make is that the genetic structure of an 
ectoproct population appears to be fundamentally like that of genetically well-known 
terrestrial animals, such as Drosophila. This conclusion is, perhaps, surprising, 
because ectoprocts are colonial, sessile, and for the most part, hermaphroditic. The 
peculiarities of the ectoproct way of life might be expected to result in unusual 
sexual and genetic mechanisms. Such mechanisms were not found. 


Breeding structure of local populations 


Local populations of Schizoporella unicornis at Green Pond and MBL dock, 
and Bugula stolonifera at Eel Pond are composed of several hundred colonies. 
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Ficure 5. Plan view and depth distribution of genotypes of “leucine” aminopeptidase 
(Locus Lap-3) in colonies of Schisoporella unicornis occurring on “islands” of pilings at the 
MBL dock. In plan view, samples are summed over all depth intervals. Genotypes appear to be 
randomly distributed, taking into account the small sample size (3-15 colonies) on each set 
of pilings. In depth, samples are summed over all areas, and by depth interval (0.0-0.6 m, 
MN = 22 00 lsm n= l4 18-30 in, 1 


Conspicuous aggregations of like genotypes are usually not found, indicating that 
larvae do not habitually settle in the immediate vicinity of the parental colony. 
Approximately random distribution of genotypes is the rule. Apparently the sessile 
mode of adult existence is not a barrier to dispersal of ectoprocts over distances of 
square meters or tens of square meters. 

Genotype frequencies at several loci in these local populations are largely in 
accord with Hardy-Weinberg equilibrium frequencies. We infer that individuals 
of each colony freely interbreed with those of other colonies in the local 
population. Since sperm and eggs are produced by the colonies of these species 
at the same time, a self-incompatibility mechanism may be present to prevent 
intra-colony fertilization. Nevertheless, some inbreeding may exist in ectoproct 
populations, but has escaped detection. As an example, the sample of 50 colonies 
of Schizoporella unicornis from the MBL dock is too small to detect a coefficient 
of inbreeding (F) of 0.25 at the 0.5 level of significance of the chi-square test. 

The close agreement with Hardy-Weinberg genotype frequencies suggests 
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that alleles at the loci sampled are selectively neutral. However, Lewontin and 
Cockerham (1958) have shown that statistical agreement with Hlardy-Weinberg 
ralues can be an insensitive indicator of selection in small population samples. It 1s 
distinctly possible that selection pressures of 10-20% are operating against particu- 
lar genotypes. There is tentative evidence for this viewpoint since genotypes of 
the Lap-3 locus in S. unicornis depart from random distribution along environ- 
mental gradients. At Green Pond a positive bias of heterozygotes appears to occur on 
the side of the floating dock closer to a source of oi] and gasoline pollution (Fig. 4+). 
At the MBL dock, the Lap-3°'/°* homozygotes become proportionately more 
abundant with increasing depth (Fig. 5). 
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Figure 6. Map showing location of collecting localities from Green Pond, MBL dock, and 
Vineyard Sound dredging, together with gene frequencies for Lap-3 alleles. Note that the 
coastal samples have nearly identical gene frequencies whereas the Vineyard Sound frequency is 
distinctly different. 
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The only locus to show statistically significant departure irom Hardy-Weinberg 
equilibrium -values (7 = 0.01 jis the 2-3 os iv Eeen oprane oi 5: 
stolonifera. There ts a marked defhciency of the heterozygote class. Since 
genotypes at the other sampled loci conform to equilibrium frequencies, self- 
fertilization or close inbreeding are not responsible. Inbreeding would lead to a 
deficiency of heterozygotes at all loci. Positive assortive fertilization could account 
for the deficit of heterozygotes, by which gametes carrying either the E-3° or 
E-3-°4 allele preferentially fertilize eggs bearing the like allele. Alternatively, 
fertilization may be random but heterozygote mortality is greater. Either explana- 
tion requires linkage equilibria between E-3 and the other loci sampled. This 
problem cannot be resolved until more is known of gamete behavior and the 
relative fitness of E-3 genotypes. 


Genetic variability of local populations 


The genome of an ectoproct probably contatns thousands of loci. Our data 
concerning genetic variability in ectoprocts are based on 8 loci in S. uwnicornis and 
11 loci in B. stolonifera, which are small samples. Other limitations of electro- 
phoretic methods in the determinations of genetic variability are discussed in Shaw 
(1965), Lewontin and Hubby (1966), and O'Brien and MacIntyre (1969). The 
most important biases recognized by these authors which are applicable to our 
research are: 


(1) A genetic variant that does not alter the mobility of its protein product 
will not be detected. 

(2) Alleles occurring at low frequencies in populations may not be sampled 
unless large numbers of individuals are analyzed. 

(3) Different mutations of the same gene may cause an identical alteration in 
mobility of a polypeptide. Thus two genetic variants may be recorded as a single 
electrophoretic variant. 

(4) Histochemical methods are known for a relatively- small number of 
proteins, most of which are enzymes. The amount of polymorphism in genes 
whose products are other classes of protein is unknown. 

Most biases inherent in electrophoresis methodology tend to underestimate the 
amount of genetic variability. Limitations (1), (2), and (3) are of this type. The 
limitations in methodology are not likely to cause recent estimates of genetic 
variability, including those reported here, to he drastically revised downward. 

Our analysis of genetic variability in ectoprocts yields estimates consistent with 
those recently made for man and several species of Drosophila (summarized by 
O’Brien and MacIntyre, 1969). The percentage of polymorphic loci in two 
ectoproct species is 25.0 and 54.5%, with a mean of 42.1%. If we add unpublished 
data for 8 loci each of two other species of Bugula, the mean value becomes 48.5%. 
The range of estimates of loci in Homo sapiens and Drosophila is 30.0-54.0%., 
with a mean of 42.0%, disregarding a single estimate of 100% based on only 3 loci. 
The average proportion of the genome heterozygous per individual (calculated 
after Lewontin and Hubby, 1966) is 0.07 in S. unicornis and 0.21 in B. stolonifera 
compared with a range of 0.12-0.33 obtained for man and Drosophila. | 

If these estimates of the amount of polymorphism are typical of the terrestrial 
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and marine taxa they represent, then little seems to distinguish the colonial, sessile, 
and hermaphroditic marine ectoprocts from solitary, free living, unisexual terrestrial 
animals in the genetic variability maintained in populations. 


Consequences of spatial separation 


The marine realm is frequently stated to be more homogenous and stable than 
the terrestrial (David, 1963) with the notable exception of the highly-variable inter- 
tidal environment (Knox, 1963). Some scientists have regarded the apparent 
scarcity of barriers to migration in the sea as evidence of sympatric speciation 
(Kohn, 1960; Weiser, 1960; and Day, 1963). Ernst Mayr, by no means an 
advocate of sympatric speciation, considered that free gene flow has led to wide- 
spread panmictic conditions in marine animals as opposed to local genetic differen- 
tiation (Mayr, 1963). 

Two contrasting models of marine ectoproct populations may be proposed. Ac- 
cording to the one just cited, extensive gene flow among populations of marine 
ectoprocts would militate against the formation of sem1-isolated local populations. 
Species would be cosmopolitan aggregations of individuals and major shifts in gene 
frequencies would only occur over extensive areas differing considerably in 
important environmental parameters such as temperature and salinity. 

The other model calls for a marine environment beset (from an ectoproct point 
of view) with barriers to gene flow.  Micro-environmental changes in substrate 
relief, texture, and composition, current strength and direction, quality and quan- 
tity of nutrients, and salinity and temperature might serve to disrupt genetic 
communication between local populations. This model would permit divergence in 
gene frequencies due either to responses to local selection pressures or, in small 
populations, to genetic drift, or a combination of selection and drift. Small, semi- 
isolated populations are commonly found in terrestrial animals such as the snail 
Cepaca nemoralis (Cain and Sheppard, 1954), species of Drosophila (Wallace, 
1968 ), and certain human populations (Cavalli-Sforza, 1959). 

Our evidence’ is insufficient to fully evaluate these models, but it is more 
consistent with a model of small, semi-isolated populations than with the concept 
of widespread panminia in the sea. An approximate 24 per cent disparity exists 
in the frequency of alleles of the Lap-3 locus in S. unicornis between the Vineyard 
Sound population and the near-shore Green Pond and MBL dock populations. 
Whether the allele frequencies change clinally or abruptly through the intervening 
populations is not yet known. In either case the 35 km interval serves as a partial 
barrier to gene flow or else stringent differential selection maintains the difference 
in the face of gene flow. Either the existence of effective barriers to gene flow or 
strong local selection pressures would disrupt panmixia among populations of 
marine species. 

More work needs to be done to assess the potentialities of dispersal in ectoprocts. 
The larvae of S. unicornis are planktonic and therefore subject to rapid current 
dispersal. However, the larvae are brooded for several davs and are released 
into the water column for a free life of only hours. They have no digestive tract 
and are nourished by a yolk supply. Estimates of current pattern and strength 
(Sumner, Osburn, Cole and Davis, 1913; Coast and Geodetic Survey tidal current 
charts, 1949) show it requires more than one day for water at the MIBL dock to 
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he transported to the Vineyard Sound location assuming a direct and rapid rate of 
transport. Thus, biological and physical limitations to rapid dispersal appear to be 
operative. Such gene flow as does occur between these populations must be via 
intervening populations. 


All of the work reported here was conducted at the Marine Biological Labora- 
tory, Woods Hole, the use of whose facilities we gratefully acknowledge. Technical 
assistance was given by Mr. 1. W. Kaufmann and'Mrs. E. Gooch. Supported in 
part by NSF erant GB=7325 te Mona [eats Scho. 


SUMMARY 


Proteins of marine species of the Phylum Ectoprocta can be easily separated by 
zone electrophoresis on polyacrylamide gels. Many bands revealed by electro- 
phoresis represent the product of individual genetic loci from which gene and 
genotype frequencies were calculated. Results are reported for 19 loci which are 
responsible for the formation of esterases, malate dehydrogenase and “lencine” 
aminopeptidase in Schisoporella unicornis and Bugula stolonifera. 

Observed genotype frequencies correspond closely with those predicted by 
Hardy-Weinberg equilibria indicating that these marine ectoprocts are dominantly 
outbreeding and that panmictic local populations cover an area of square meters. 
In addition, 25.0-54.5% of the diagnosed loci are polymorphic, depending on the 
species. Forty-two % of all loci identified are polymorphic. 

The most important implication of these findings is that the population structure 
and amount of genetic variability in ectoprocts (and by implication in other marine 
organisms) is fundamentally similar to that of terrestrial organisms. If these 
tentative conclusions are substantiated by subsequent work, the processes of specia- 
tion in marine organisms may differ little from those proposed in the better docu- 
mented research on terrestrial animals. 


PERATURE CITED 


ALLEN, S. L., 1961. Genetic control of the esterases in the protozoan Tetrahymena pyriformis. 
Ann. New York Acad. Sci., 94: 753-773. 

Ar-Rusuo1, A. H., 1958. The polygenic basis of sex-ratio in Tigriopus. Proe. 10th Int. 
Congr. Genet., 2: 8-9. 

Avucustinsson, K-B., 1961. Multiple forms of esterase in vertebrate blood plasma. Ann. New 
York Acad. Sci., 94: 844-860. 

Baricozzi1, C., 1960. Genetic systems in sessile and semisessile organisms, pp. 559-569. In: 
A. A. Buzzati-Traverso, Ed. Perspectives in Marine Biology. University of 
California Press, Berkeley. 

Batractta, B., 1958. Balanced polymorphism in Tisbe reticulata, a marine copepod. Evolution, 
12: 358-364. 

Battrac.ia, B., 1965. Advances and problems of ecological genetics in marine animals. Proe. 
llth Int. Congr. Genct., 2: 451-463. 

BECKMAN, L., AnD F. M. JouNson, 1964a. Esterase variations in Drosophila melanogaster. 
Heriditas, 51: 212-220. 

BECKMAN, L., AND F. M. Jonnson, 1964b. Genetic control of aminopeptidases in Drosophila 
melanogaster. Hceriditas, 51: 221-230. 

Bocguet, C., 1951. Recherches sur Tisbe (-Jdyea) reticulata, n. sp. Arch. Zool. Exp. Gén., 
87: 335-416. 


CAA 
Cal 


CENETIC OF BC LOPROCTS 1 


BocgueEt, C., 1953. Recherches sur le polymorphisme naturel des Jaera martina (Fabr.) 
(Isopodes Asellotes). Arch. Zool. Exp. Gén., 90: 187-450. 

BULLIVANT, J. S., 1967. Release of sperm by Bryozoa. Ophelha, 4: 139-182. 

Burns, J. M., ann F. M. Jounson, 1967. Esterase polymorphism in natural populations of 
a sulfur butterfly Colias curytheme. Science, 156: 93-96. 

BuzzatTi-Traverso, A. A., 1960. Perspectives in marine biology, pp. 613-621. Jn: A. A. 
Buzzati-Traverso, Ed., Perspectives in Wlarme Biology. University of California 
Press, Berkeley. 

CAVALLI-SForzA, L. L., 1959. Some data on the genetic structure of human populations. Proc. 
10th Int. Congr. Genet., 1: 389-407. 

Carn, A. J., AND P. M. Suepparp, 1954. The theory of adaptive polymorphism. Amer. 
Natur., 88: 321-326. 

Coast AnD GEODETIC Survey, 1949. Tidal current charts: Narragansett Bay to Nantucket 
Sound, second edition. Serial No. 628. U. S. Dept. Commerce, 12 charts. 

Davin, P. M., 1963. Some aspects of speciation in the Chaetognatha, pp. 129-143. Ju: J. P. 
Harding and N. Tebble, Eds., Speciation m the Sea. The Systematics Association, 
London. 

Day, J. H., 1963. The complexity of the biotic environment, pp. 31-49. In: J. P. Harding 
and N. Tebble, Eds., Speciation in the Sca. The Systematics Association, London. 

FALCONER, D. S., 1960. Introduction to Quantitative Genetics. Ronald Press Co., New 
kork 305 pp. 

FoS B. ann [. Gygen, 1954. Studies on Serpulid Pomotoceros triquetcr L. IT. The 
colour pattern of the branchial crown and its inheritance. Nytt Jlag. Zool., 2: 85-90. 

Gorm A. M., 1967. Sulfonylation with sulfonyl halides, pp. 706-711. In: C. H. W. Hirs, Ed, 
Methods in Ensymology, Volume II, Enzyme structure. Academic Press, New York. 

iene eet... «1963. Protein differences in Drosophila. I. Drosophila melanogaster. 
Genctics, 48: 871-879. 

Hussy, J. L.. ann R. C. Lewontin, 1966. A molecular approach to the study of genic hetero- 
zygosity in natural populations. I. The number of alleles at different loci in 
Drosophila pseudoobscura. Genetics, 54: 577-594. 

Hyman, L. H, 1959. The Invertebrates, F'olume 5. McGraw-Hill Book Co, New York, 
783 pp. 

Jacosson, K. B., 1968. Alcohol dehydrogenase of Drosophila: interconversion of isoenzymes. 
Science, 159: 324-325. 

Jonna FOX, C. G. Kaxarr, R. H. RICHARDSON, M. R. WHEELER AND WS. STONE, 1966, 
An analysis of polymorphisms among isozyme loci in dark and light Drosophila 
ananassae strains from American and Western Samoa. Proc. Nat. Acad. Sci. U.S.A. 
56: 119-125. 

Knox, G. A., 1963. Problems of speciation in intertidal animals with special reference to 
New Zealand shores, pp. 7-29. Im: J. P. Harding and N. Tebble, Eds., Speciation 
in the Sea. The Systematics Association, London. 

Konv, A. J., 1960. Problems of speciation in marine invertebrates, pp. 571-588. In: A. A. 
Buzzati-Traverso, Ed., Perspectives in Marine Biology. University of California 
Press, Berkeley. 

Kojima, K., anp K. M. Yarsroucu, 1967. ‘Frequency dependent selection at the Estcrase-6 
locus in Drosophila melanogaster. Proc. Nat. Acad. Sci. U. S. A., 57: 645-649. 

LEWONTIN, R. C., AnD C. C. CocKERHAM, 1958. The goodness-of-fit test for detecting natural 
selection in random mating populations. Evolution, 13: 561-564. 

LewonTIn, R. C, ano J. L. Hussy, 1966. A molecular approach to the study of genic 


heterozygosity in natural populations. II. Amount of variation and degree of 
heterozygosity in natural populations of Drosophila pseudoobscura. Genetics, 54: 
595-609. 


MAKINO, S. 1951. an Atlas of the Chromosome Numbers in Animals. Towa State College 
Press, Ames, Iowa, 240 pp. 

MANWELL, C., 1966. Sea cucumber sibling species: Polypeptide chain types and oxygen 
equilibrium of hemoglobin. Science, 152: 1393-1396. 

MANWELL, C., AND C. M. A. BAxer, 1963. A sibling species of sea cucumber discovered by 
starch gel electrophoresis. Comp. Biochem. Physiol., 10: 39-53. 


156 JAMES EYGOOCH ANXN DTFEONAS TAOL sChOrr 


MANWELL, C.. AND C. M. A. BAKER, 1968. Genetic variation of isocitrate, malate and 
6-phosphogluconate dehydrogenases in snails of the genus Cepaeua—introgressive hybrid- 
ization, polymorphism and pollution? Comp. Biochem. Physiol., 26: 195-209. 

MANWELL, C., C. M. A. Baxer, P. A. ASHTON AND ETD S. Corner, 1967. Biochemical 
differences between Calanus finmarchicus and C. helgolandicus. J. Mar. Biol. Ass. 
U. K., 47: 145-169. 

Maturo, F. J. S., Jr., 1966. Bryozoa of the southeast coast of the United States: Bugulidae 
and Beaniidae (Cheilostomata: Anasca). Bull. Mar. Sci., 16: 556-583. 

Mayr, E., 1963. Animal Species and Evolution. Belknap Press of Harvard University 
Press, Cambridge, Massachusetts, 797 pp. 

MILKMAN, R., 1967. Genetic and developmental studies on Botryllus schlosseri. Biol. Bull., 
132: 229-243. 

Moore, R. O., AND C. A. ViLLEE, 1963. Multiple molecular forms of malic dehydrogenase in 
echinoderm embryos. Comp. Biochem. Physiol., 9: 81-94. 

O'Brien, S. J., AND R. J. MacIntyre, 1969. An analysis of gene-enzyme variability in natural 
populations of Drosophila melanogaster and D. simulans. Amer. Natur., 103: 97-113. 

Ocita, Z, AND T. Kasar, 1965. Genetico-biochemical analysis of specific esterases in Musca 
domestica. Jap. J. Genet., 40: 173-184. 

Ray, D. L., 1960. Some marine invertebrates useful for genetic research, pp. 497-512. In: 
A. A. Buzzati-Traverso, Ed., Perspectives in Marine Biology. University of California 
Press, Berkeley. 

RyLanp, J. S., 1960. The British species of Bugula (Polyzoa). Proc. Zool. Soc. London, 
134: 65-105. 

RyLanpb, J. S., 1965. Catalogue of Mam Warie Fouling Organisms, Vol. 2: Polysoa. 
Organization for Economic Co-operation and Development, Paris, 83 pp. 

SABADDIN, A., 1959. Analisi genetica del policromatismo di Botryllus schlosseri (Pallas) 
Savigny (Ascidiacea). Boll. Zool., 26: 221-243. 

SCANDALIOS, J. G., 1969. Genetic control of multiple molecular forms of enzymes in plants: 
A review. Biochem. Genet., 3: 37-79. 

ScHopr, T. J. M., ann F. T. MANnHEIM, 1967. Chemical composition of Ectoprocta (Bryozoa). 
J. PaléOnéol., 31: 1197-1225. 

Scuopr, T. J. M., and F. T. Manuetm, 1968. Chemical composition of Ectoprocta (Bryozoa) : 
Addenda. J. Paleontol., 42: 858. 

Sexton, E. W. ANd A. R. Crark, 1936. Summary of the work on the amphipod Gummarus 
chevreuxi Sexton carried out at the Plymouth Laboratory (1912-1936). J. Mar. Biol. 
Ass. Us KM2 TSA. 

SHaw, C. R. 1965. Electrophoretic variation in enzymes. Science, 149: 936-943. 

SILEN, L., 1966. On the fertilization problem in the gymnolaematous Bryozoa. Ophelia, 13: 
113-140. 

STAIGER, H., 1954. Der Chromosomendimorphismus beim Prosobranchier Purpura lapillus in 
Beziehung zur Ökologie der Art. Chromosoma, 6: 419-478. 

STONE, W. S., M. W. WHEELER, F. M, JoHNSon anp K. Nojima, 1968. Genetic variation in 
natural island populations of Drosophila nasuta and Drosophila ananassae subgroups. 
Proc. Nat. Acad. Sci. U. S. A., 59: 102-109. 

SUMNER, F. B., R. C. Ospurn, L. J. Core anp B. M. Davis, 1913. -A Biological Survey of 
the Waters of Woods Hole and Vicinity. Aull. Bur. Fish. U. S., 31 (2 pts): 1-861. 

Tuornser, E. J., I. T. OLtver ann P. B. Scurr, 1968. Comparative electrophoretic patterns 
of dehydrogenases in different species. Comp. Biochem. Physiol., 25: 973-987. 

Wacrace, B., 1968. Topics in Population Genetics. W. W. Norton and Co., Inc., New York, 
481 pp. 

Wieser, W., 1960. Problems of species formation in benthic microfauna of the deep-sea, pp. 
513-518. Jn: A. A. Buzzati-Traverso, Ed., Perspectives m Mlarine Biology. Univer- 
sity of California Press, Berkeley. 

Wricut, T. R. F., 1963. The genetics of an esterase in Drosophila melanogaster. Genetics, 
18: 787-801. 

YARBROUGH, K., AND K. Koyrma, 1967. The mode of selection at the polymorphic Esterase-6 
locus in cage populations of Drosophila melanogaster. Genetics, 57: 677-686. 





